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Atom-dimer scattering below the three-body break-up threshold is studied for a system of three 
identical bosons. The atom-dimer scattering length and the energy of the most weakly-bound three- 
body state are shown to be strongly correlated. An appropriate rescaling of the observables reveals 
the subtlety of the correlation, and serves to identify universal trends in the unitary limit of divergent 
two-body scattering length. The correlation provides a new quantitative measure of the degree of 
universality in three-body systems with short-ranged interactions, as well as a consistency check of 
effective field theories and other theoretical models. 
PACS numbers: 



Among the most striking results of recent experimen- 
tation with supercooled atomic gases is the demonstra- 
tion that trap loss rates are extraordinarily sensitive to 
few-body interactions within a trapped many-body sys- 
tem. The experiments of many groups [lj show signa- 
tures of few-body correlations within a trapped ensem- 
ble of Bose alkali vapors at nano-Kelvin temperatures. 
This discovery has stimulated a large number of theoret- 
ical and experimental investigations, particularly in the 
unitarity limit of a divergent two-body scattering length, 
where the system has no natural length scale beyond that 
of the trapping potential. 

It is widely anticipated that in the vicinity of unitar- 
ity, supercooled Bose gases display universal collective 
properties. However, criteria for the onset of universal- 
ity and measures of the degree of universality have not 
yet emerged. The primary purpose of this Letter is to 
provide such a measure through a thorough investiga- 
tion of atom-dimer scattering for a wide variety of two- 
body short-ranged potentials and the zero-range interac- 
tion model. 

Weakly-bound few-body systems with relatively large 
two-body scattering lengths have long been studied in 
nuclear physics. A salient example is a result published 
by Phillips in 1968 which puzzled nuclear theorists for 
more than a decade. Phillips compared results of calcu- 
lations of the neutron-deuteron scattering length and the 
energy of the triton bound state made with different two- 
body potentials. He found that, unlike two-body scat- 
tering in which the energy of the last bound state scales 
as the inverse squared scattering length (E ~ — 1/a 2 ), 
the bound state energy of the triton (n-n-p) appears 
approximately proportional to the scattering length of 
the neutron-deuteron collision. The data illustrating this 
correlation is shown in Figure [1] 

Many subsequent works Q confirmed that a strong 
linear correlation between these two observables exists. 
Two decades later a simple and elegant explanation of 
the Phillips line was given by Efimov and Tkachenko 
[3] , who simply noted that the apparent linearity derives 
from the fact that the two-body potentials used sample 
only a small portion of the space of scattering parameters; 
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FIG. 1. The Phillips line, showing the unexpected linear 
correlation between the triton bound state energy and the 
neutron-deuteron scattering length (data from ref. Q). 



that is, they yield similar values for the neutron-deuteron 
scattering lengths. 

While the origin of Phillips' observation is now well un- 
derstood, modern experiments with ultracold gases are 
able to sample a much wider range of scattering param- 
eters by magnetic field tuning through a Fano-Feshbach 
resonance. This suggests that the correlation between 
two- and three-body parameters can be investigated in 
far greater detail than in earlier nuclear physics stud- 
ies. Theoretical and experimental works related to uni- 
versality in ultracold Bose gases have been concentrated 
on three-body recombination in the close vicinity of the 
three-body threshold. In this work we discuss the proper- 
ties of the three-body system at the two-body threshold 
and below. 

In this work we introduce a new parametrization of 
the relationship between the three-body (atom-dimer) 
scattering length and the energy of the last three-body 
bound state for the specific case of three identical Bosons. 
This relationship is referred to below as the modified 
Phillips line. In contrast to the original Phillips line, 
this relationship is found to be linear over a large range 
of interaction parameters and more directly reflects the 
well known threshold law for single-channel scattering, 
namely E ps It also provides a simple test for uni- 

versality of three-body systems. 
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FIG. 2. a) Phillips line for three 4 He atoms: the Hez bind- 
ing energy as a function of the He — He2 scattering length 
based on four commonly used two-body potentials. Note the 
approximately linear relationship (dashed line), b) Modified 
Phillips line for three 4 He atoms for the same set of potentials. 
Note an improved linearity. 



As pointed out by Efimov and Tkachenko [4(, for a 
weakly-bound three-body state the following approxi- 
mate correlation should hold 
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where E 2 is the dimer binding energy and m\ 2 is the re- 
duced mass of the particle-dimer system. A more trans- 
parent representation of the strength of the correlation 
is obtained by rewriting Eq. [T] in terms of the scaled 
dimensionless variables 
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The variable a can be thought of as a dimensionless scat- 
tering length, and uj characterizes the three-body binding 
energy. If the three-body state nearest to the threshold 
is deeply bound, w is small; large values of ui indicate the 
existence of a weakly-bound three-body state. 

Using a recently developed three-body code, we have 
strenuously tested the well studied case of three bosonic 
4 He atoms @, and so we will use the helium trimer as a 
first illustration of our improved parameterization. The 
traditional Phillips line for helium trimer states calcu- 
lated with four commonly used two-body potentials is 
plotted in Figure [2^,, where £3 is the energy of the He 3 
threshold bound state (relative to the three-body break- 
up threshold) , and 03 is the scattering length for a helium 
atom and a bound helium dimer. Figure ^jp shows the 
modified Phillips line obtained using the suggested scal- 
ing with the dimer binding energy. Note that the scaling 
results in an improved fit. 

Equation [1] derives from the fact that atom-dimer scat- 
tering is dominated by the pole of the t-matrix corre- 
sponding to the near-threshold state of the trimer. Gen- 
erally speaking, we should expect the linear relation of 
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FIG. 3. (Color online) A modified Phillips line plotted with 
rescaled parameters for a wide variety of two-body parameters 
and several alternative two-body potentials. Red crosses mark 
the results for realistic He-He potentials (see Fig. [2j 



the rescaled parameters a and oj to hold providing that 
no other poles of the i-matrix are relevant in the en- 
ergy range of interest. We have tested this assertion for 
systems with three indistinguishable bosons by studying 
the correlation for a variety of two-body potentials with 
widely adjusted two-body scattering lengths and binding 
energies. 

Figure[3]illustrates the correlation for six different fam- 
ilies of two-body potentials. The potentials include the 
TTY potential with an artificial coupling constant, the 
family of Bargmann potentials [f| with fixed small ef- 
fective range and a scattering length varying from 2 to 
250 atomic units, the Bargmann potential with an ef- 
fective range simulating the He-He interaction, the fam- 
ily of Bargmann potentials with varying asymptotic nor- 
malizing constant, and the MTV potential ("symmetric 
model" for the triton) with a varying coupling constant. 
The near-linearity over a wide range of u is apparent 
in the figure, though deviations can be seen, especially 
for small values of the 3— body scattering length, mag- 
nified in the inset of Figure |3l where the complexity of 
the correlation is revealed. An interesting aspect of the 
plot, which we had not anticipated, is the emergence of 
a universal behavior as the 3— body scattering length ap- 
proaches zero. 

To better understand the complex correlation revealed 
in Figure [3j we show in Figure |4] the same correlation 
plot using a single two-body potential (the Bargmann 
potential with the effective range ro = 1 a.u.), but with 
large variations of the two-body scattering length. We 
have also shown the three-body data generated by solv- 
ing the regularized Skornyakov/Ter-Martirosyan equa- 
tions (STM) 0,i] corresponding to a zero-range interac- 
tion model. 

We can identify two distinctive regimes: a universal 
regime and a strong coupling regime. In the strong cou- 
pling regime (characterized by relatively small two-body 
scattering length) the shape of the modified Phillips line 
depends on the potential (inset Fig. [3] ) and, typically, 
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forms an elongated half-loop above the universal curve. 
First, as the two-body scattering length increases, the 
points on the Phillips line move right along the correla- 
tion plot reaching a local maximum, and then turn left 
as they approach the universal regime. 

In the universal regime, as the two-body scattering 
length increases, the points along the Phillips line move 
left along the universal curve. The shape of the univer- 
sal correlation plot stems from the interplay of poles of 
the three-body t-matrix corresponding to the formation 
of near-threshold bound and virtual states. 

We can distinguish three characteristic parts of the uni- 
versal curve. The linear part at the right end corresponds 
to large positive particle-dimer scattering length and a 
very shallow three-body bound state. The diverging part 
at the left end of the correlation plot corresponds to large 
negative particle-dimer scattering length and indicates 
the formation of a near-threshold virtual state. With 
increasing two-body scattering length this virtual state 
turns into a bound state, and the corresponding point 
on the modified Phillips line jumps to the far end of the 
positive linear part of the universal curve. There is also 
a transitional nonlinear regime, when the scaled particle- 
dimer scattering length a is small, and both poles con- 
tribute to the shape of the universal curve. 

The most transparent interpretation for the shape of 
the universal part of the correlation plot can be obtained 
from a simple analysis of the spectrum of the Faddeev 
operator 

K{E) = G 2 {E)V{P + + P~) , 

where G 2 {E) = (H + V - E)- 1 is the Green's func- 
tion for the three-body cluster Hamiltonian, Ho is the 
Hamiltonian for three free particles, V is the two-body 
pairwise potential and P ± are Jacobi coordinate trans- 
formation operators. The equation for the component of 
the scattering wave function then reads 

[1 + K{E)]$ = -V(P+ + P-) X0 , 

where xo stands for the atom-dimer plane wave. The 
component of the bound state wave function satisfies the 
homogeneous equation 

[1 + = . 

The Faddeev operator K(E) for short-range potentials 
has a discrete spectrum with eigenvalues X n (E). For E < 
Ei the eigenvalues X n (E) are real. 

At the two-body threshold E — E 2 the Faddeev op- 
erator can be approximated by a sum of projectors on 



1 The three-body wave function can be recovered from the com- 
ponent as * = (1 + P + + P - )*. 
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FIG. 4. The modified Phillips line for the Bargmann potential 
and the zero-range model. For the Bargmann potential the 
effective range is held fixed at 1 a.u., while the scattering 
length varies from 2 to 460 a.u. . The values of the two-body 
scattering length 02 are shown along the curve. 

the states corresponding to eigenvalues A + and A~ clos- 
est to —1, and the particle-dimer scattering length can 
be expressed as 



d l + \+(E 2 ) l + \-(E 2 ) ■ v ; 

There are two possible situations: 1) one of the eigen- 
values - X + (E 2 ) or \~(E 2 ) - is very close to the critical 
value A = — 1 and the corresponding term gives the ma- 
jor contribution to the particle-dimer scattering length; 
2) both of the terms contribute substantially. The first 
case is responsible for the "linear" and "singular" parts 
of the a(uj) universal correlation plot. The second case 
corresponds to the intermediate regime of small scaled 
particle-dimer scattering length a. 

This interpretation is illustrated in Fig. [31 where we 
show the particle-dimer scattering length (top) and the 
three-body operator spectrum (bottom) as a function 
of the two-body scattering length. The poles in the 
three-body scattering length correspond to the eigenvalue 
of the three-body operator crossing the critical value 
\(E 2 ) = —1, where a new three-body bound state is 
formed. 

The universal part of the modified Phillips line is re- 
produced by the simple formula 

a = - — + uj + a Q (4) 

cj u;o 

with a\ = 5.5, ujq = 0.419 and ctQ — 4. The first term 
here is responsible for the description of the large nega- 
tive scattering length regime, the other two terms fit the 
large positive scattering length regime. Equation 2] can 
be used in practical calculations to estimate the trimer 
binding energy from the scattering length. A simpler 
empirical fit a = a; + |(1 — 1/w 2 ) is suitable for positive 
three-body scattering length. 

Careful analysis of the modified Phillips line reveals 
that in the universal, potential-independent regime the 
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few minutes, depending on the desired numerical accu- 
racy and physical parameters of the system). 

The code is available from the authors by request and 
will be available online in the near future. 
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FIG. 5. Atom-dimer scattering length and the spectrum of 
the kernel of the Faddeev (a) or STM (b) equation at the 
two-body threshold as a function of the two-body scattering 
length, a) Bargmann potential with ro = 1 a.u., b) zero-range 
interaction model. 



scattering length for the atom-dimer collision approaches 
zero when E3/E2 « 2.54. The other special point on the 
plot corresponds to divergent negative scattering length, 
which corresponds to the universal ratio E3/E2 ~ 6.7. 

The result demonstrates a novel form of universality 
in weakly bound systems. Simple in physical nature, the 
modified Phillips line provides an important test for nu- 
merical and theoretical analysis of 3-body systems. It 
affords the opportunity to check bound-state and scat- 
tering results for consistency and classifies 3-body sys- 
tems according to distinct dynamical regimes. It also 
provides an opportunity to check estimates of 3-body 
bound states and atom-dimer scattering lengths for inter- 
nal consistency; such as those obtained with regularized 
zero-range potential models Q . Correlation plots similar 
to the modified Phillips line for bosons can be constructed 
for three-body systems with nontrivial spin-isospin struc- 
ture and nonidentical particles. Here we shall only men- 
tion that the data shown in Figure[T]is consistent with the 
universal part of the modified Phillips line constructed on 
the base of the zero-range interaction model Q . 

All the calculations presented here have been per- 
formed using an original code for solving Faddeev equa- 
lol - 12j • The code is currently capable of three-body 



tions 

bound state and scattering calculations below the first 
excitation threshold. It provides users a simple interface 
for preparing configuration files and for constructing op- 
timized grids automatically [lOj ■ The computational ker- 
nel solves the system of Faddeev equations for bound or 
scattering states. It uses a numerically effective compu- 
tational scheme which - in combination with an option 
of using multiple CPUs - makes it possible to calculate 
scattering amplitudes quickly (from a few seconds to a 
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